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Abstract 

We report the first observation of the decay J/ip — > 37. The signal has a statistical significance 
of 6a and corresponds to a branching fraction of B(J/ip — » 37) = (1.2 ± 0.3 ± 0.2) x 10 -5 , in which 
the errors are statistical and systematic, respectively. The measurement uses ip(2S) — ► J/ip 
events acquired with the CLEO-c detector operating at the CESR e + e~ collider. 

PACS numbers: 13.20.Gd, 12.38.Qk 
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Ortho-positronium (oPs), the 3 Si e + e~ bound state, decays to 37 almost exclusively and 
has long been a fertile ground for precision QED tests [l|. The analog to oPs— > 37 for 
quantum chromodynamics (QCD), three-photon vector quarkonium decay, has not yet been 
observed. The rate of three-photon J / if) decays acts as a probe of the strong interaction 0] , 
most effectively when expressed in relation to J /if) — > jgg, J /if) — > 3g, or J /if) — > due 
to similarities at the parton level. Hence, measurements of £> 37 , B igg , B 3g , and £>« relative to 
one another (where Bx = B(J/if) — > X)) provide crucial experimental grounding for QCD 
predictions @, 0, 0] . 

In this Letter we report the first observation of J/if> — > 37. Rate measurements for other 
rare or forbidden all-photon decays, J /if) — > 77, 47, 57, and •yr] c with n c — > 77, are also 
described. Previous searches for uj and J /if) decay to 37 have yielded branching fraction 
upper limits of 1.9xl0~ 4 and 5.5xl0~ 5 , respectively As with oPs, C-parity symmetry 
suppresses vector quarkonia decays to an even number of photons, and two-photon decays are 
forbidden by Yang's theorem 0). Ref. 0] reports the limit £> 77 <2.2xl0~ 5 at 90% confidence 
level (C.L). Five-photon decays are suppressed by an additional factor of (at least) ~ a 2 ; 
c.f., £(oPs -> 57) ~ 2 x 1(T 6 Q . 

Ignoring QCD corrections altogether, Ref. 0] predicts £>3 7 /£>« ~ a/14, B^/B igg ~ 
(a/a s ) 2 /3 and B^/B^g ~ (a/a s ) 3 . Using the precisely measured Bu p] in the first pre- 
diction implies £?3 7 ~ 3xl0~ 5 . The latter two suffer the uncertainty of what value of a s 
to employ at the charmed quark mass scale 0. Assuming a s (m 2 )=0.3 and inserting the 
result from a recent CLEO measurement [s| {B igg ~ 0.09 and B 3g ~ 0.66) into the latter two 
predictions gives B 3l w (0.9-1.6)xl0 -5 . The first-order perturbative QCD corrections j3] to 
these estimates are large, so these predictions should only be considered as approximate. 

Events were acquired at the CESR e + e~ collider with the CLEO detector [lOj , mostly in 
the CLEO-c configuration (95%) with the balance from CLEO III. The dataset corresponds 
to 27xl0 6 produced ip(2S) mesons and (9.59±0.07)xl0 6 ip(2S) -> n + n-J/ip decays 
Event selection requires the tracking system to find exactly two oppositely charged particles, 
corresponding to the ir + 7T~ recoiling from the J /if), and that the calorimeter have at least 
2, 3, 4, 5, and 3 photon showers for the J /if) — > 77, 37, 47, 57, and 7^ c (— ► 77) samples, 
respectively. Photon candidates must have energy exceeding 36 MeV and, with respect 
to any shower associated with one of the charged pions, either be located (a) more than 
30 cm away, or (b) between 15 cm and 30 cm from it and have a photon-like lateral shower 
profile. We require that photon candidates not be located near the projection of either 
pion's trajectory into the calorimeter nor be aligned with the initial momentum of either 
pion within 100 mrad. 

A two-step kinematic fit first constrains the beam spot and the two charged pion candi- 
dates to a common vertex, and then the vertexed 7r + 7r~ and the most energetic n photon 
candidates to the if)(2S) mass [E| and initial three-momentum, including the effect of the 
~3 mrad crossing angle between the e + and e~ beams. Tight quality restrictions are ap- 
plied to the vertex (x 2 /d.o.f.<3) and four- momentum (% 2 /d.o.f.<3) fits. The mass recoiling 
against the 7r + 7r~ must lie inside a window around the J /if) mass, M(ir + 7T~ — recoil) =3087- 
3107 MeV. Non- J /if) backgrounds are estimated by keeping a separate tally of events with 
M(7T + n~— recoil) inside 2980-3080 MeV or 3114-3214 MeV, ranges which together are ten 
times wider than the signal window. 

Events with any of the photon pairs in the mass windows 0.10-0.16 GeV, 0.50-0.60 GeV, 
or 0.90-1.00 GeV are rejected to eliminate contributions from decays with 7r°'s, rfs, or 77 "s, 
the dominant sources of photons in J / if) decays. For the 37 selection only, we require all pho- 
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ton pair masses be less than 2.8 GeV to eliminate potential contamination from rj c — > 77. 
This requirement effectively restricts the smallest energy photon to have energy exceed- 
ing 200 MeV. For the 47 and 57 samples only, the smallest shower energy must be above 
120 MeV, and all lateral shower profiles must be photon-like. This last restriction on shower 
shape avoids feed-up from J/ip — > 777C), vO^ll events with one or more photon conver- 
sions between the tracking chambers and the calorimeter: in such cases the two showers 
from the conversion e + and e~ overlap one another, thereby distorting both of their lateral 
profiles. For the 777c channel only, we restrict the search region to large M(77)i g and small 
^(77)sm, which are, respectively, the largest and smallest of the three two-photon mass 
combinations in the event. The signal region is chosen this way so as to keep backgrounds 
small. Specifically, the signal box is defined, in units of GeV, by 0.16<M(77) sm <0.48, 
2.985<M(77) lg +0.0935M(7 7 ) sm <3.040. 

Signal and background decay modes are modeled with Monte Carlo (MC) samples that 
were generated using the EvtGen event generator 12], fed through a GEANT-based 13| 



detector simulation, and then exposed to event selection criteria. For J ftp — > nj signal 
decays, final state photon momenta are distributed according to phase-space. For J/tp — > 37, 
the lowest order matrix element for ortho-positronium [14j is used as an alternate; compared 
to phase space, it modestly magnifies the configurations that are two-body-like and those 
with three nearly-equal-energy photons (at the expense of topologies lying between these two 
extremes). For the process J/if> — > 7^ c , an n c mass and width of 2979.8 MeV and 27 MeV, 
respectively, are used (both are close to the PDG values to generate a Breit-Wigner 
77-mass distribution; alternate widths from 23-36 MeV and different lineshapes [I5J are 
explored as systematic variations. 

Distributions in M(77) sm vs. M(77)i g and M{n + it~ — recoil) for the J /ip — > 37 and J/tj) — > 
7^(77) samples are shown for data, signal MC samples, and likely background decays in 
Figs. [1] and [2J respectively. 

In all modes, non-J/ip backgrounds are small and are subtracted statistically using 
M(-k + tt~— recoil) sidebands in the data. We determine the backgrounds from J/ip decays 
with an exhaustive study of Monte Carlo samples. Decays with J/ip — > 7/j (where fj sig- 
nifies any of the many isoscalar mesons in the mass range from 600-2500 MeV), followed by 
fj ~ > 77 pose a negligible threat for any of the target modes because the product branching 
fractions are extremely small (e.g., ~2xl0 -8 for J/ip — * 7/ 2 (1270), / 2 (1270) — > 77). The 
predominant source of backgrounds to the 37 sample is the 77r°7r° final state. This type of 
event can survive the selection by having both 7r° decay axes nearly parallel to their lines of 
flight, such that one photon of each pair has very low energy in the laboratory frame, and 
is therefore nearly irrelevant to conservation of four-momentum. An analysis by BES [l6[ 
found that the largest sources of J/tj) — > 77r°7r° are from J/tj) — > jfj decays, specifically 
through / 2 (1270) and / (2050), followed in importance by / (1710), / (1500), and a number 
of much smaller contributions from nearby resonances. However, not all relevant product 
branching fractions for J/ip — > jfj, fj — > 7r°7r° have been measured, those that are measured 
have large uncertainties, and interference effects among overlapping fj may not be small. A 
method to normalize 77r°7r° other than using measured branching fractions is employed to 
reduce systematic uncertainty. The x 2 /d.o.f. distribution for 77r°7r° decays has a character- 
istic shape, nearly independent of 7i°7i° mass, as shown in Fig. [3j the region x 2 /d.o.f.=5-20, 
where almost no signal is present, is used to establish the level of J/ip — > 77r°7r°. Figure H] 
shows the x 2 /d-o.f. distribution from data, MC signal and MC background and the small 
contribution from non-J/^ decays obtained from the M(tt + 7t~ — recoil) sidebands. The 142 
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data events with x 2 /d.o.f.=5-20 contain J /if) — > 37 signal (3.4 events), non-J/if) background 
(3.2), and, using known branching fractions, J /if) — > u)i], 77 — > 77 (1.7), J /if) — > 777, 77 — > 77 
(1.2), J/if) -> 777, 77 -> 3vr° (1.2), J/V -> 777', 77' -> 7CJ, u) -» tt°7 (0.6), and J/V> -> 7vr° 
(0.2). The remainder (130.5 events) serves to normalize the 77r°7r° background component, 
which has a relative 8% statistical uncertainty. 

As a cross check on the 37 background normalization, we perform a maximum likelihood 
fit to data in the entire J /if) — > 37 x 2 /d.o.f .=0-30 region with the combination of shapes from 
MC of 77r°7r° and 37 signal with floating normalizations for each, and a fixed J / ^-sidebands 
contribution from data, scaled by a factor of 0.1. Using this method with different sources 
of the 77r°7r° taken one at a time as 100% of the background results in an average signal 
size of 23.3 events (with variation from 22.8 to 24.1), which is 0.9 events smaller than our 
nominal technique. Based on these numbers we assign a systematic error of 0.9 events, or 
~5% relative, for signal extraction and background estimation for J/if> — > 37. 

The x 2 /d.o.f. fit just described is repeated with the 37 signal shape weight fixed to zero. 
The likelihood difference with respect to the nominal fit provides a measure of the statistical 
significance of the signal. This significance varies from 5.9a to 6.6a when using any one of 
the backgrounds 7/ 2 (1270), 7/ (1500), 7/ (1710), 7/ (2020), 77r°7r° (phase space) as the 
sole contributor to the background shape. 

MC studies indicate the following primary sources of backgrounds for the other modes: 
for the 27 sample, J /if) — > 771" (3.3 events) and 77/, 77 — > 77 (2.7); for the 47 sample, 
J /if) — > 777, 77 — > 77 (0.9) or 77 — > 3n° (0.8), 777', 77' — > 77 (0.3) or 77' — > ju), u — > 7r°7 
(0.9) or 77' — > 7r°7r°77, 77 — > 77 (0.3); for the 57 sample, J /if) — *> 777, 77 — > 3tt° (0.2) and 777', 
77' — > 7r°7r°77, 77 — > 77 (0.3); for 777c J/if> ^ 777, 77 — > 77 (0.3), 777', 77' — > 77 (0.2), and our 
newly found signal, <//■?/> — >■ 37 (0.3). 

Numerical results appear in Table [B Net yield uncertainties and upper limits on event 
counts include the effects of statistical fluctuations in signal and background estimates. 
Signal efficiencies range from ~2% (57) to ~22% (37), and J/if> — > 37 is the only mode 
with a clear signal: 37 events observed on a background of 12.8. Statistics dominate the 
overall uncertainties for all decay modes. The J /if) — > 37 efficiencies for pure phase-space 
and the oPs matrix element are equal to within (0.2±0.1)%; nevertheless, a 15% systematic 
error is assigned to allow for different behavior in the much heavier J /if) system. For 777,,, 
uncertainties in the lineshape, background, and T(r] c ) dominate the systematic error. 

Using the recently determined B(J/if) —>■ 7^ c )= (1.98±0.09±0.30)% [l5j], the 77 c — > 77 
branching fraction can be calculated as B(i] c -> 77 ) = (0.6±S;i±0.1) xlO" 4 , or <3xl0" 4 at 
90% C.L. This value is consistent with the PDG fit value of (2.7±0.9) x 10~ 4 at the level of 
1.3cr, although making a meaninful comparison is difficult because the PDG number depends 
indirectly upon previous, considerably smaller values for B(J/if> —> 777c). 

In conclusion, we have investigated decays J /if) —> 727 with n = 2,3,4,5, where the 
photons are produced in direct decay, not through an intermediate resonance. For n=3, a 
signal of 6a significance is found with branching fraction ^3 7 =(1.2±0.3±0.2) xlO -5 . This 
value lies between the zeroth order predictions [ij for B^/B igg and B 3l /B 3g and is consistent 
with both, but is a factor of ~2.5 below that of B&y/Bu- This measurement represents the 
first observation of a three-photon meson decay. No signal is seen for n—2, 4, or 5, and 
upper limits are set on the branching fractions, each of which is the most precise or only 
measurement. We also measure B(J/ip —> 7r/ c ) x £>(t7 c — > 77)= (1.21^^0.3) xlO -6 or an 
upper limit of <6xl0 -6 at 90% C.L., both consistent with other determinations [5j. 
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TABLE I: Results for the five J/ip — > wy decay modes, showing the raw number of signal can- 
didate events, estimated background levels, statistical significance of each signal, the net event 
yield and its 90% C.L. upper limit (UL), the signal efficiency, different sources of systematic error 
and their quadrature sum, expressed in percent of the central value (37, jTj c ) or of the UL (oth- 
ers), the branching fraction B(J/ip — > X) [first error is statistical, second is systematic], and the 
corresponding 90% C.L. upper limit UL, including effects of systematic errors. 
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FIG. 1: Top four plots: in 37 data (lower left) and MC events for different J/tp decays (top row 
and lower right), the largest vs. the smallest two photon mass combination per event. In the MC 
plots, darker shading of each bin signifies higher event density than lighter shading; in the data 
plot, each dot represents an event. The solid lines demarcate regions excluded from the J ftp — > 37 
selection. Bottom plot: distribution of M(ir + ir~ — recoil) for the data events (points with error bars) 
overlaid with the J ftp — > 37 signal MC prediction (dotted line histogram) and MC background plus 
signal (solid line histogram) normalized to the data population. The arrows indicate the region of 
accepted recoil mass. 
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FIG. 2: As in Fig. [H except zoomed in on the rj c region, and the overlaid parallelogram indicates 
the J /ip — > 77/ c signal region. The M(tt + tt^— recoil) distribution is for 77/ c candidates. 
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FIG. 4: The distribution of x 2 /d.o.f . for J/ip — > 37, in the top plot showing data (points with error 
bars) overlaid with the sum (dotted line histogram) of three components: non- J/ip background from 
scaled data sidebands (shaded histogram) and MC predictions for signal (solid) and J/ip — ► 771-° 7r° 
background (dashed). The bottom plot shows the same distribution, but with the MC and non- 
J / ip background subtracted from the data. The arrows indicate the values for signal selection and 
background normalization. 
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